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Abstract 

The NGC 1333 IRAS 4A protobinary was observed in the 1.3 cm and 6.9 mm continuum and the ammonia and 
SiO lines, with an angular resolution of about 0.4 arcseconds. The continuum maps show the circumstellar structures 
of the twoprotostars, Al and A2. The Al system is brighter and more massive than the A2 system. The ratio of mass, 
including dense gas and protostar, is about 6. The properties of the circumstellar disks and outflows suggest that Al 
may be younger than A2. The deflected part of the northeastern jet of A2 is bright in the SiO line, and the distance 
between the brightest peak and deflection point suggests that the enhancement of SiO takes about 100 yr after the 
collision with a dense core. The ammonia maps show a small structure that seems to be a part of the obstructing 
core. The outflow properties were studied by comparing interferometric maps of SiO, ammonia, formaldehyde, and 
HCN lines. Their overall structures agree well, suggesting that these species are excited by the same mechanism. 
However, the intensity distributions show that SiO is chemically unique. SiO may be directly linked to the primary 
jet while the other species may be tracing the entrained ambient gas. 

Key words: ISM: individual (NGC 1333 IRAS 4A) — ISM: jets and outflows — ISM: structure — stars: 
binaries: general — stars: formation 
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1. INTRODUCTION 

NGC 1333 IRAS 4A2 is one of the best-studied low-mass 
protostars. It has a circumstellar disk viewed nearly edge-on 
and drives a well-collimated and rotating outflow/jet (Choi 
2005; Choi et al. 2007, 2011). Choi et al. (2010) found that 
the NH3 (2, 2) and (3, 3) lines seem to selectively trace the ac- 
cretion disk while the millimeter/centimeter continuum seems 
to trace both disk and envelope. The rotation kinematics of the 
disk suggests that the mass of the central protostar is ^^0.08 
Mq and that the collapse age is ^50,000 yr (Choi et al. 
2010). The bipolar jet flows in the northeast-southwest direc- 
tion (Blake et al. 1995; Hodapp & Ladd 1995; Lefloch et al. 
1998). The jet may be launched from a small region on the 
disk, or outflow foot-ring, with a radius of ^2 AU (Choi et 
al. 2011). The northeastern jet shows a sharp bend of flow di- 
rection that may be caused by a collision with a dense core in 
the ambient cloud (Choi 2005; Back et al. 2009). IRAS 4A2 
also exhibits other star formation activities such as H2O maser 
emission (Park & Choi 2007; Marvel et al. 2008). 

IRAS 4A2 belongs to a binary system (IRAS 4A) in the 
Perseus star-forming region at a distance of 235 pc from the 
Sun (Lay et al. 1995; Looney et al. 2000; Hirota et al. 2008). 
The spectral energy distribution of IRAS 4A suggests that the 
constituent binary members are Class protostars (Sandell et 
al. 1991; Enoch et al. 2009). The luminosity of IRAS 4A sug- 
gests that the protostars are growing at a rate typical of Sun-like 
stars through the accretion of gas from the molecular envelope 
(Jennings et al. 1987; Choi et al. 2010). 



IRAS 4A1 is brighter in the millimeter/centimeter contin- 
uum than its sibling protostar A2 (Looney et al. 2000; Reipurth 
et al. 2002; J0rgensen et al. 2007), but it is weaker in the NH3 
emission (Choi et al. 2007). IRAS 4A1 seems to drive an out- 
flow to the south, but its counter flow, if any, has not been de- 
tected clearly (Choi 2005). While Al must be in an evolution- 
ary stage similar to that of A2, the details are less clear. 

In this paper, we present the results of our observations of the 
NGC 1333 IRAS 4 A region in the 1.3 cm and 6.9 mm contin- 
uum, the NH3 lines, and the SiO line with an angular resolution 
higher than those of the previous studies (Choi 2005; Choi et 
al. 2007). We describe our observations in Section 2. In Section 
3, we report the results and discuss the star-forming activities 
in the IRAS 4A region. A summary is given in Section 4. 

2. OBSERVATIONS AND DATA 

2. 1. K-band Observations 

The NGC 1333 IRAS 4 region was observed using the 
Very Large Array (VLA) of the National Radio Astronomy 
Observatory in the NH3 (2, 2) and (3, 3) lines and in the A 
= 1.3 cm continuum. Observations in the C-array configuration 
were presented in Choi et al. (2007). Details of the observations 
in the B-array configuration and results for the NH3 lines were 
given in Choi et al. (2010). The continuum was observed for 
15 minutes at the beginning and for another 15 minutes at the 
end of each observing track, and the NH3 lines were observed 
in the midsection of the track. 

For the 1.3 cm continuum, the observations were made in 
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Fig. 1. Maps of the continuum emission toward the NGC 1333 IRAS 4A region. The contour levels are 1, 2, 4, 8, 16, 32, and 64 X 0.05 mjy beam"^. 
Dashed contours are for negative levels, (a) Map of the A = 1.3 cm continuum made using the B-array data only. The rms noise is 0.016 mJy beam~^. 
Shown in the bottom right-hand comer is the synthesized beam: FWHM = 0''40 X 0"33 with PA = 64°. (b) Map of the A = 6.9 mm continuum. A 
relatively bright substructure is labeled as A3. The rms noise is 0.017 mJy beam^^. The synthesized beam has FWHM = 0."58 X Of 52 with PA = -1°. 
Plus signs: the 3.6 cm continuum sources (Reipurth et al. 2002). 



the standard K-band continuum mode (22.5 GHz). Maps were 
made using a CLEAN algorithm. With a natural weighting, the 
B-array continuum data produced a synthesized beam of 0"40 
X 0."33 in full-width at half-maximum (FWHM). 

2.2. Q-band Observations 

The NGC 1333 IRAS 4 region was observed using the VLA 
in the SiO u = J = 1 ^ line (43.423858 GHz) and in the 
A = 6.9 mm continuum. Twenty-five antennas were used in the 
C-array configuration on 2004 March 1. The continuum was 
observed for 20 minutes at the beginning and for 10 minutes at 
the end of each observing track, and the SiO line was observed 
in the midsection of the track. For the SiO line, the spectral 
windows were set to have 64 channels with a channel width 
of 0.098 MHz, giving a velocity resolution of 0.67 km s^^. 
For the 6.9 mm continuum, the observations were made in the 
standard Q-band continuum mode (43.3 GHz). 

The phase tracking center was (a, 5) = (03''29"U0?993, 
31°13'45."21) in J2000.0. Note that the field of view of the C- 
array observations corresponds to the first field of the D-array 
mosaic observations (Choi 2005), which covers the northeast- 
ern part of the IRAS 4A area including the redshifted outflow 
of IRAS 4A2. The southwestern field was not covered in the 
C-array observations, and most part of the blueshifted outflows 
are outside the field of view. 

The nearby quasar 0336H-323 (PKS 0333H-321) was observed 
to determine the phase and to obtain the bandpass response. 
The flux was calibrated by observing the quasar 0713H-438 
(QSO B0710H-439) and by setting its flux density to 0.20 Jy 
(VLA Calibrator Flux Density Database''). To check the flux 
scale, the quasar 0319H-415 (3C 84) was observed in the same 
track. Comparison of the amplitude gave a flux density of 4.28 



See http://aips2.nrao.edu/vWcalflux.html. 



Jy for 0319H-415, which agrees with the value in the VLA 
Calibrator Flux Density Database. The bootstrapped flux den- 
sity of 0336H-323 was 1 .26 Jy. To avoid the degradation of sen- 
sitivity owing to pointing errors, pointing was referenced by 
observing the calibrators at the X band (A = 3.6 cm). This ref- 
erenced pointing was performed about once an hour and just 
before observing the flux calibrator With a natural weighting, 
the SiO data produced a synthesized beam of FWHM = 0''55 
X 0''52, and the continuum data produced a beam of FWHM 
= 0"58 X 0."52. 

3. RESULTS AND DISCUSSION 

3.1. Continuum Emission Maps 

Figure 1 shows the continuum emission maps, and source 
parameters are listed in Table 1. The L3 cm continuum map 
shows the compact structures of IRAS 4A1 and A2 as well as 
a few weak features. IRAS 4A1 is marginally resolved, and 
an elliptical Gaussian fit gives a deconvolved size of FWHM 
= 0"24 X 0;'16 with a position angle (PA) of 52°. IRAS 4A2 
is unresolved. The 6.9 mm continuum map shows IRAS 4A1, 
A2, and an extended structure surrounding the binary system. 
IRAS 4A2 does not stand out because the extended structure 
has a comparable brightness. Some of the emission peaks in 
the extended structure are located away from Al or A2 with 
angular distances comparable to the A1-A2 binary separation, 
which suggests that there is a common envelope around the 
binary system. 

The measurement of total flux is a nontrivial issue, especially 
for the 6.9 mm map, because the source structure is compli- 
cated. The total flux of Al is not very sensitive to the choice of 
integration area (a box in this case) because the main peak of 
Al is much brighter than the other features in the area. The to- 
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Table 1. NGC 1333 IRAS 4A Continuum Source Parameters 





Peak Position" 


1.3 cm Flux Density'' 


6.9 mm Flux Density'' 


Source 


aj2000.0 <5j2000.0 


Peak Total 


Peak 


Total 


Al 

A2 


03 29 10.53 3113 31.0 
03 29 10.42 3113 32.1 


1.554 ±0.016 2.47 ±0.07 
0.225 ± 0.016 0.37 ± 0.07 


5.57 ± 0.02 
0.39 ± 0.02 


9.51 ±0.12 
1.46 ±0.08 


A 




2.66 ±0.18 




11.5 ±0.3 



° Source positions are from the 1.3 cm continuum map. Units of right ascension are hours, minutes, and 
seconds, and units of declination are degrees, arcminutes, and arcseconds. 

'' Flux densities are in mJy beam~^ or in mJy, corrected for the primary beam response. The box sizes used 
for measuring the total flux are the same as those listed in Table 2 of Looney et al. (2000): 2"9 x 2"2 for 
Al, l."8 X l."6 for A2, and 5. "4 x 6. "2 for A. 



tal flux of A2 at 6.9 mm, however, is sensitive to the choice of 
box because of the extended emission around A2. Since com- 
parisons with the 2.7 mm flux densities are important (see be- 
low), we measured the total flux densities (listed in Table 1) 
in the boxes used by Looney et al. (2000), for consistency in 
measuring spectral index. However, the Al and A2 boxes have 
an overlapping region, and the A box is larger than the union 
of Al and A2 boxes. Therefore, the flux density of A is not a 
simple sum of those of Al and A2. In the calculation of mass 
given below, we mostly rely on the 1.3 cm and 2.7 mm total 
flux densities of the whole IRAS 4A region, and the ambiguity 
in the 6.9 mm flux density of A2 is not critical. 

Figure 2 shows the spectral energy distributions (SEDs) in 
the centimeter-millimeter wavelength region. SEDs can be de- 
scribed using a power-law form, F oc i/" , where F is the flux 
density, ly is the frequency, and a is the spectral index. Dust 
emission from a disk/envelope system usually has a large spec- 
tral index (a ^ 2), and free-free emission from a thermal radio 
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Fig. 2. Spectral energy distributions. Filled squares: IRAS 4A (includ- 
ing Al, A2, and common envelope). Open squares: IRAS 4A1. Open 
circles: IRAS 4A2. Flux densities are from Rodriguez et al. (1999), 
Reipurth et al. (2002), this work, and Looney et al. (2000). The flux 
uncertainties are smaller than the size of markers. 



jet usually has a small spectral index (a ^ 1) (Reynolds 1986; 
Anglada et al. 1998). The SED of the IRAS 4A system shows 
that the slope is steep at shorter wavelengths and shallow at 
longer wavelengths. The spectral index is large (a J^ 1.9) at 
wavelengths shorter than 3.5 cm, and small (a w 0.6) in the 
6.2-3.5 cm section. This change of a around 3.5 cm suggests 
that the continuum emission at 1 .3 cm and shorter wavelengths 
mostly comes from dust and that the 6.2 cm continuum may 
come from free electrons. Even at 3.5 cm, a significant fraction 
of the flux may come from dust. 

The mass of molecular gas can be estimated from the SED 
of the dust continuum using the mass emissivity given by 
Beckwith & Sargent (1991), 
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vn 



/3 



cm g 



(1) 



where i^g = 1200 GHz, and /3 is the opacity index. For the 1.3 
cm - 2.7 mm section, the opacity index of IRAS 4 A is /? w 1 .3. 
Assuming optically thin emission, the mass can be estimated 
by 

A/ = ^f^, (2) 

where D is the distance to the source. By is the Planck func- 
tion, and Td is the dust temperature. 

Assuming D = 235 pc and T^ = 33 K (Hirota et al. 2008; 
Jennings et al. 1987), the mass of IRAS 4A is M = 2.6 ± 0.8 
A/0. The mass estimate is sensitive to the value of /?. For ex- 
ample, Looney et al. (2000) derived a mass of 1.3 Mq (scaled 
to the distance of 235 pc) assuming (3=1 and using the 2.7 
mm flux density only. Therefore, the uncertainty in mass in- 
troduced by the assumptions on /3 may be about a factor of 2. 
This mass is based on the flux densities in a box of 5 "4 x 6 "2 
around IRAS 4A1, which may include the circumstellar disk 
and protostellar envelope around Al, those around A2, and the 
common envelope. For comparison, the mass of the whole pro- 
tostellar envelope of IRAS 4A derived from single-dish obser- 
vations is r^l Mq (Enoch et al. 2009, scaled to 235 pc), and the 
mass of the IRAS 4A2 protostar derived from the kinematics of 
die circumsteflar disk is -0.08 Mq (Choi et al. 2010). 

The mass of molecular gas around each protostar can be de- 
rived using the ratio of flux densities, assuming that the dust 
properties are nearly uniform in the whole IRAS 4A system. 
IRAS 4A1 and A2 are most clearly separated in the 1.3 cm 
continuum map, and the 1.3 cm flux densities give 2.3 ± 0.7 



M. Choi et al. 



Mq for Al and 0.3 ± 0.1 A/© for A2. Then the total mass 
of the A2 system (protostar and dense molecular gas) is '^0.4 
Mq. Therefore, regardless of the mass of the Al protostar, the 
Al system is much more massive than the A2 system. 

The extended structure in the 6.9 mm continuum map 
(Figure lb) is clumpy and shows several substructures. The 
overall morphology of the extended structure is similar to what 
was seen in the 2.7 mm continuum (see Figure 12d of Looney 
et al. 2000). One of the substructures is relatively bright and 
can be seen in both the 6.9 mm and 2.7 mm maps. This sub- 
structure, labeled as A3 in Figure lb, has a peak position at 
0''7 west and 0''7 south with respect to Al. The nature of the 
substructures is not clear. Some of them may be dense clumps 
in the envelope (that may eventually fall into one of the pro- 
tostars), some may be prestellar cores, and some might even 
harbor small compact objects. 

3.2. NH3 Line Maps of the IRAS 4A1 Region 

Figure 3 shows high-resolution (B-array) NH3 maps of the 
region around IRAS 4A1, and Figure 4 shows the spectra. (See 
Choi et al. 2007 for NH3 maps in a lower resolution.) The 
source structure and spectra of Al is somewhat different from 
those of A2 (see Choi et al. 2010 for high-resolution maps of 
the A2 region). In the case of A2, the NH3 lines seem to trace 
the circumstellar disk since the source shape is elongated in 
the direction perpendicular to the outflow and the spectra has 
a profile consistent with a rotating disk (Choi et al. 2010). By 
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Fig. 3. Maps of the NH3 lines toward IRAS 4A1 from the VLA 
B-an'ay data. Contribution from the continuum emission was subtracted 
out. (a) Map of the (2, 2) line core, Vlsr = (5-8, 7.6) km s^^. (b) Maps 
of the blueshifted, (4.5, 5.8) km s"!, and redshifted, (7.6, 8.9) km s"\ 
parts of the (2, 2) line, (c) Map of the (3, 3) line core, (5.8, 7.6) km s^^. 
(d) Maps of the blueshifted, (4.6, 5.8) km s" \ and redshifted, (7.6, 8.8) 
km s~^, parts of the (3, 3) line. Contour levels are 3, 4, and 5 times the 
rms noise [0.35 mjy beam^^ for (a) and (c), and 0.43 mJy beam^'^ 
for (6) and (d)] . Shown in the bottom right-hand corner is the restoring 
beam: FWHM = 0!'3. Cross: peak position of the 1.3 cm continuum. 
Arrow: direction of the southern outflow (toward outflow peak 16). 
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Fig. 4. Spectra of the NH3 lines toward IRAS 4A1, at the peak posi- 
tion of the 1.3 cm continuum emission. Vertical dotted line: systemic 
velocity of IRAS 4A (Vlsr = 6.7 km s-i; Choi 2001). 

contrast, the NH3 maps of Al show a more complicated mor- 
phology, and the spectra of Al show relatively weak line cores 
and wider line wings than those of A2. The total integrated 
NH3 flux of Al is smaller than that of A2 by a factor of 1 .7. 

The (2, 2) line core map (Figure 3a) shows a structure asso- 
ciated with the continuum source. This structure (central com- 
ponent) is elongated in the east-west direction, which may be 
the circumstellar disk, considering that the blueshifted outflow 
of Al is directed toward the south (Choi 2005). The (3, 3) fine 
maps (Figures 3c and 3d) seem to show the central compo- 
nent and a secondary structure about 0."4 south (southern com- 
ponent). The southern component is blueshifted and may be 
related with the southern outflow. The peak positions of the 
blueshifted and redshifted line wings of the central component 
(Figure 3d) could not be separated with the 0''3 beam, and the 
rotation of the Al disk, if any, may be slower than that of the 
A2 disk. Therefore, the rotation curve could not be analyzed 
with our data. 

In summary, the circumstellar structure of Al traced by the 
NH3 lines seems to be a mixture of the accretion disk and the 
southern outflow. The Al disk is weaker than the A2 disk in 
the NH3 lines. The rotation of Al disk seems to be slower than 
that of A2, suggesting that the Al protostar is less massive than 
the A2 protostar 

5.5. Evolutionary Status of the IRAS 4A System 

The IRAS 4A system is one of the youngest and best-studied 
protostar binary systems. The evolutionary status of A2 is es- 
pecially well-studied thanks to the rotation kinematics of the 
circumstellar disk (Choi et al. 2010). An interesting issue is the 
evolutionary status of Al in comparison with A2. 

Previously, Choi et al. (2007) examined the IRAS 4A system 
with NH3 and continuum images of ^1" resolution. To under- 
stand the difference in the NH3-to-continuum flux ratios they 
explored two possibilities: (1) the protostars are exactly coeval, 
or (2) A2 is a protostar, and Al may be a transitionary object 
on the verge of protostellar collapse. Because the second pos- 
sibility contradicts other observational facts, they preferred the 



NGC 1333 IRAS 4A Imaging 



first possibility. However, it implied that the dust-to-gas ratios 
of Al and A2 may be quite different, by almost an order of 
magnitude, which is not easy to understand. 

The new higher-resolution images presented in this paper 
suggest that the difference in the properties of the disks may 
be the key in understanding the NHs-to-continuum flux ratios 
and also the evolutionary status of the IRAS 4A system. Here 
we revisit the issue and explore the possibility that Al and A2 
are protostars at slightly different stages of evolution. That is, 
though Al and A2 are members of a single binary system, their 
collapse ages can be different, i.e., not exactly coeval. Then 
which one is the younger of the two? We suggest that Al is 
younger than A2 for several reasons. 

First, the accretion disk of A2 traced in the NH3 lines shows 
a clearly detectable rotation kinematics, while the Al disk does 
not. Since the disk around more massive object rotates faster, 
the A2 protostar may be more massive than the Al protostar. 
Therefore, if the accretion rates are similar, the Al protostar 
may be younger than the A2 protostar. (If the accretion rate of 
Al is higher, as the envelope of Al is more massive, then the 
age of Al may be even much smaller than that of A2.) 

Second, the A2 disk is brighter than the Al disk in the NH3 
lines. The strength of the NH3 lines is likely correlated with the 
mass of the disk, and the disk mass increases monotonically in 
the simplest models of protostellar evolution (Young & Evans 
2005). Therefore, the difference in the NH3 line strength im- 
plies that, if the accretion rates are similar, the Al disk is less 
massive and younger than the A2 disk. 

Third, the outflow of A2 is much longer than that of Al. 
The length of the northeast-southwest bipolar outflow driven 
by A2 is --130" (from the A2 protostar to the H2 knot HL3; 
see Figure 1 of Choi et al. 2006). The length of the southern 
outflow driven by Al is ^15" (from the Al protostar to the 
SiO outflow peak 16; see Figure 2a of Choi 2005). Assuming 
a proper motion of 0''064 yr^^ (Choi et al. 2006) for rough 
estimates, the time scale of the Al outflow is shorter than that 
of A2by-l,800yr 

The relative youth of Al, despite its relatively massive 
circumstellar material, suggests that the progenitor prestellar 
cores evolved at different speeds. As a simplest scenario of 
the evolution of binary prestellar cores, suppose that the par- 
ent cloud fragmented into two cores: core 1 with a mass of 
Mci that later evolves to become the Al system, and core 
2 with a mass of Mc2 to become the A2 system. If the col- 
lapse age of A2 is indeed larger than that of Al, core 2 may 
have evolved faster and formed a protostar earlier. Since the 
denser core evolves faster (that is, the evolution time scale is 
shorter for the denser core), core 2 may have had a higher den- 
sity than core 1 . This scenario can be understood with a sim- 
ple argument. At the time of the fragmentation, suppose that 
each of the two cores has a roughly spherical shape (or any 
shape similar to each other), and the boundary (the point of 
contact between the two spheres) is the point where the grav- 
itational forces of the two cores balan ce out. Then the ratio 
of radii is Rci/Rc2 ~ \/Mci/Mc2, and the density ratio 
is nci/nc2 ~ ^/\/Mci/Mc2- Though the exact density ra- 
tio may depend on details of the fragmentation process, the 
more massive core may usually have a smaller average density 



than the sibling core. (For the IRAS 4A system, if the mass ra- 
tio at the time of fragmentation is similar to the current value, 
Mci/Mc2 ~ Mai/Ma2 ~ 6 and nci/nc2 ~ 0.4.) Therefore, 
it is reasonable that, when a cloud fragments into two cores, 
the less massive one would be denser, evolve faster, and form a 
protostar earlier. Considering the amount of circumstellar ma- 
terial, however, Al may outgrow A2 in the future and eventu- 
ally become a more massive star 

If Al is indeed younger than A2, how large is the age gap? 
To answer this question, the collapse age of Al needs to be 
measured using the rotation kinematics of the disk, as what 
was done with A2, which requires imaging with a resolution 
higher than ours. Instead, here we simply calculate a charac- 
teristic time scale of the system as a sketchy expectation based 
on simple physics. The age gap may be in the order of a time 
scale that governs the dynamics of the A1-A2 system, i.e., a 
sound-crossing time. Assuming an effective sound speed of 
0.27 km s^^ (Young & Evans 2005), the projected separation 
of I'.'S (430 AU) gives a sound-crossing time of ^7,500 yr. 
Considering that the A1-A2 system in not necessarily in the 
plane of the sky, the value above is a lower limit. For example, 
if we assume that the A1-A2 system is ^30° out of the plane 
of the sky, the sound-crossing time would be ~9,000 yr, which 
is ~20% of the collapse age of A2. (See Figure 5 for the com- 
parison between the binary separation and the infall region.) 
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Fig. 5. Schematic diagram of the IRAS 4A system. Protostars Al, A2, 
and BI are labeled. Plus signs: SiO outflow peaks. See Figure 2a of 
Choi (2005). Dashed lines: outflow axes. Gray circle: infall region with 
a radius of 2,700 AU, which corresponds to the size of the spherical 
collapse wave with a sound speed of 0.27 km s^^ and a collapse age 
of 50,000 yr 
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[See http://minho.kasi. ie.kiyPublications.html for the original figure.] 
Fig. 6. Maps of the NH3 (3, 3) line toward the IRAS 4A region, made from the VLA C-array data (also see Figure 1 of Choi et al. 2007). The maps 
were convolved to have an angular resolution of FWHM = I'.'O as shown in the bottom left-hand comer. For these maps, the continuum levels were not 
subtracted from the line data, and the structure associated with IRAS 4A1, A2, and BI comes from both the continuum and line emission, (a) Map of the 
Une core. Black contours show the intensity distribution averaged over the velocity interval of Vlsr = (5.8, 7.6) km s^^. Contour levels are 3, 4, and 5 
times the rms noise (0.7 mJy beam^^). (b) Maps of the hne wings. Blue and red contours are for the velocity intervals of (-1.0, 5.8) and (7.6, 14.4) km 
s~^, respectively. Contour levels are 3, 4, 5, 6, 7, and 8 times the rms noise (0.4 mJy beam^^). Some outflow peaks are labeled (also see Figure 5). Gray 
scale: map of the SiO v = Q J = I —^ Q line (Figure 3a of Choi 2005). Plus signs: the 3.6 cm continuum sources (Reipurth et al. 2002). 




Fig. 7. Spectra of the NH3 (3, 3) line from the C-array maps with an angular resolution of 2"0 (Figure 6) toward several NH3 outflow peaks and IRAS 
4A1/2. Vertical dotted hne: systemic velocity of IRAS 4A. 



3.4. Outflows used, the maps were convolved with larger beams to show the 

„, , , r, • 1 TTi A n .< » • J J outflows clearly. Figure 5 shows a schematic diagram of the 

The molecular outflows m the IRAS 4A region were detected t^. » f. ^ « it- ^ , »t.t /o os i- 

1 xTT. /I OS 1- 1 1 1 11 /o OS 1- IRAS 4A system, and Figure 6 shows NH3 (3, 3) line maps 

in the NH3 (3, 3) line but almost undetected in the (2, 2) line. . , , , • .c o// r- -^ 1 / 

„. ... J. , „ ' with an angular resolution of 2 . Figure 7 shows the spectra at 

Since some weak emission structures of the outflows tend to „ , . . „. „ , , xtt, ,^ on T 

, , . ,, . , , . , , . ,, outflow peak positions. Figure 8 shows the NH-? (3, 3) line map 
be buried by noise when the nominal synthesized beams were 
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Fig. 8. Maps of the northeastern area of the IRAS 4A region. The con- 
tinuum levels were not subtracted from the line data. Black contours: 
map of the NH3 (3, 3) line from the VLA C-array data, averaged over 
the velocity interval of Vlsr = (6.4, 7.6) km s^^. Contour levels are 3, 
4, and 5 times the rms noise (0.5 mjy beam^^ ). Negative levels are not 
shown, to avoid clutter. Shown in the bottom left-hand corner is the syn- 
thesized beam: FWHM = 0."98 X 0."95 with PA = 47°. Red contours: 
map of the SiO t; = OJ=l— >0 line from the VLA C-array data, av- 
eraged over the velocity interval of (7.7, 25.9) km s^^. Contour levels 
are 1, 2, 3, 4, and 5 X 0.5 mJy beam^^, and the rms noise is 0.18 mJy 
beam^^. The SiO map was convolved to have an angular resolution of 
FWHM = l."0, Protostars and outflow peaks are labeled. 

with an angular resolution of 1" for the velocity range show- 
ing the outflow peak (OP) 3a clearly. Figure 8 also shows the 
SiO line map showing the northeastern jet of IRAS 4A2. For 
these maps (Figures 6 and 8), contributions from the contin- 
uum emission were not removed because the velocity extents 
of the line wings at strong peaks are comparable to the spectral 
coverage and left only a little room for line-free channels that 
can be used for determining the continuum levels. The struc- 
tures directly associated with the continuum sources in these 
maps should be ignored. 

In general, most of the NH3 peaks associated with the IRAS 
4A1/2 outflows are located on or near the SiO outflow peaks 
(Figure 6b), but the intensity distributions are quite different. In 
the SiO outflow maps, the strongest redshifted emission comes 
from OP 1 and 2, and the strongest blueshifted emission comes 
from OP 7 and 9. By contrast, in the NH3 maps, OP 3 and 4 
are the strongest redshifted peaks, and OP 15 and 16 are the 
strongest blueshifted peaks. 

OP 17 is unusually strong in the NH3 line, but the nature 
of this source is not clear (Lefloch et al. 1998; Choi 2005). In 
addition to the IRAS 4A outflows, the NH3 map also shows the 
outflow associated with IRAS 4BI (the redshifted emission at 



~4" north of BI in Figure 6b). 

3.5. Deflection of the Northeastern Outflow 

One of the most interesting properties of the IRAS 4A2 
outflow is the directional variability of the northeastern (red- 
shifted) jet. The change of flow direction, or a sharp bend, near 
OP 3 is a good example of jet deflection caused by a collision 
with a dense core in the ambient cloud (Choi 2005). Figure 8 
shows the detailed structure of the SiO jet, and Figure 9 shows 
the spectra. The undeflected upstream part (OP 3-5) is weak in 
the SiO line. The deflected downstream part (OP 1-3) is rela- 
tively strong in the SiO line. OP 2 in the low-resolution map 
(Choi 2005) is now resolved into two peaks (OP 2a and 2b) 
in the new map. Though these two peaks are spatially close to 
each other (within 2."5 or 600 AU), their velocity structures are 
very different: OP 2a has double peaks (one at 17 and the other 
at 21 km s^^), and OP 2b has a single peak (at 18 km s~^). This 
sudden change of velocity structure suggests that the deflected 
jet is highly turbulent. 

One of the motivations of the new high-resolution observa- 
tions was imaging the detailed structure of the jet at the jet-core 
impact point that is supposed to be near OP 3. However, the 
SiO emission of OP 3 is weak and does not show a compact 
structure, which may mean that the SiO emission does not get 
bright immediately after the impact but takes some time to de- 
velop. A likely explanation is that Si atoms in dust grains are 
injected into the shocked gas at the impact point, combine with 
oxygen in the gas, and form SiO in a time scale of 100-1,000 yr 
(Schilke et al. 1997; Gusdorf et al. 2008; Guillet et al. 2009). 
The SiO molecules then combine again with oxygen and be- 
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Fig. 9. Spectra of the SiO line toward SiO outflow peaks 1-3 
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[See http://minho.kasi. re. ki"/Publications. html for the original figure.] 
Fig. 10. Maps of the molecular outflows in the IRAS 4A region. For these maps, continuum levels were subtracted out. (a) Maps of the H2CO 141 GHz 
line (Choi et al. 2004). Blue and red contours show the intensity distiibution averaged over the velocity intervals of Vlsr = (-1.0, 5.8) and (7.7, 14.4) 
km s~^, respectively. The lowest contour level and the contour interval are 70 mJy beam^^, and the rms noise is 22 mjy beam~^. Negative levels are 
not shown, to avoid clutter (b) Maps of the HCN J = 1 — >■ line (Choi 2001). See Figures 1 and 2 of Choi (2001) for the velocity intervals. Note that 
the blueshifted and redshifted hne wings cannot be clearly separated because of the blending of hyperfine components. The lowest contour level and the 
contour interval are 60 mJy beam^^, and the rms noise is 20 mJy beam^^. Gray scale: map of the SiO line (Figure 3a of Choi 2005). Plus signs: the 3.6 
cm continuum sources (Reipurth et al. 2002). 



come Si02. Therefore, the strength of the SiO emission would 
peak with a certain time lag after the jet-core impact. 

In the deflected jet of IRAS 4A2, the SiO emission is bright- 
est at OP 2a, ^7" from the sharp bend. There is no measure- 
ment of the proper motion in this part of the outflow. If we 
use the proper motion measured in the southwestern jet (0."064 
yr^^; Choi et al. 2006) as a very rough guide, the time scale of 
outflow propagation from the impact point to OP 2a is ^100 
yr. The argument above assumes that the jet is a steady flow, 
but images of the IRAS 4A2 jet show that the brightness dis- 
tribution is far from uniform. The brightness distribution of the 
deflected jet would also be subject to the non-steady structure 
of the impacting jet. Considering the average separation be- 
tween outflow peaks, the uncertainty in the time scale of SiO 
peak emission would be ^40%. 

The NH3 maps (Figures 6 and 8) show an interesting struc- 
ture around OP 3. While the emission peaks directly associ- 
ated with the outflow (OP 3b and 3c) show wide redshifted 
line wings, those located slightly away from the SiO jet (OP 3a 
and 3d) show narrow lines near the ambient velocity (Figure 
7). OP 3a and 3d are probably parts of the dense core that 
is obstructing the northeastern jet and causing the deflection. 
NH3 molecules at OP 3a is excited either by the shock propa- 
gated from the impact point or by a stream of jet that penetrated 
deep into the obstructing core. Examinations of the centimeter- 
millimeter continuum maps show that there is no detectable 
continuum source at or near the position of OP 3a (Rodriguez 



et al. 1999; Reipurth et al. 2002; Choi et al. 2004, 2007; this 
work). 

3.6. Outflow Chemistry 

As mentioned above, there is an interesting difference be- 
tween the intensity distributions of SiO and NH3. Comparisons 
with other outflow tracers would be helpful in understand- 
ing the intensity distributions. Figure 10 shows the H2CO and 
HCN maps of the IRAS 4A outflows (Choi 2001; Choi et al. 
2004). The outflow structure in SiO is most well-collimated 
and straight, and the outflow images in the other tracers agree 
with the SiO image in the overall structure. Outflow peaks 
in the four tracers are usually located close to each other. 
Especially, OP 1, 3, 4, and 15 are detectable in all the four out- 
flow maps, and the peak positions of each OP in these maps are 
clustered within ^2". (The average separation between adja- 
cent peaks is ~6".) This clustering suggests that a single phys- 
ical agent (e.g., internal shock in the primary jet) is responsible 
for the excitation of all the four outflow tracers. Therefore, the 
outflowing gas components traced by these tracers are physi- 
cally closely connected to each other and flow with a certain 
degree of coherence. 

Though the outflow tracers share the source of excitation, 
SiO stands out when the intensity distributions are compared. 
In the northern part of the IRAS 4A region, the SiO map shows 
that the deflected (downstream) part of the northeastern jet (OP 
1-3) is brighter than the undeflected (upstream) part (OP 3-5). 
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Fig. 11. Correlation diagrams of the SiO, NH3, H2CO, and HCN line intensities for outflow peaks (see the maps in Figures 6b. 10a, and 10b). The 
integrated intensities are corrected for the primary beam response. Uncertainties are marked for outflow peaks 1, 3, 4, and 15. Written in the top/bottom 
right-hand corners are the linear congelation coefficients. 



The maps of the other tracers show that the upstream part is 
brighter In the southern part of the region, the SiO map shows 
that the southwestern jet (blueshifted jet of IRAS 4A2: OP 7- 
12) is brighter than the southern outflow (blueshifted outflow 
of IRAS 4A1: OP 13-16). The maps of the other tracers show 
that the southern outflow (especially OP 15 and 16) is brighter 
The intensity distribution pattern of the NH3, H2CO, and HCN 
maps is similar to that of the CO J = 2 -^ 1 outflow map (see 
Figure 1 of Girart et al. 1999). 

To show this trend quantitatively, correlation diagrams of the 
outflow peaks were made. The intensity of the SiO line tends 
to be anti-correlated with those of the other lines (Figure 11, 
top panels). The linear correlation coefficients are negative: 
r « —0.3. If only the outflow peaks detected in all the four lines 
are considered (marked with uncertainty bars in Figure 11), the 
anti-correlation is even stronger: r w —0.9. By contrast, the 
intensities are positively correlated among NH3, H2CO, and 
HCN: r « 0.9 (Figure 11, bottom panels). This difference sug- 
gests that the outflow component traced by SiO is fundamen- 
tally different from that of the other tracers. 

Blake et al. (1995) arrived at a similar conclusion based on 
the spectral profiles of molecular lines obtained from single- 
dish observations, which show that the SiO emission is kine- 
matically displaced from the bulk cloud velocity. The image 
and spectra of the SiO line suggests that the gas traced by the 
SiO emission is either the primary jet itself or an outflow com- 



ponent closely linked to the primary jet, as suggested by the ob- 
servations of other outflows such as L1448 and HH 212 (Dutrey 
et al. 1997; Codefla et al. 2007). 

The dichotomy between SiO and the other outflow tracers 
can be explained in two different ways. First, the differences in 
shock strength may trigger different chemical processes. While 
the sputtering of Si from dust grains requires shock velocities 
higher than ^25 km s~^, species such as H2CO do not sur- 
vive at such high velocities (Blake et al. 1995; Schilke et al. 
1997; Arce et al. 2007). In this scenario, the intensity ratio 
of SiO to the other species can be an indicator of the shock 
strength. Second, the line intensities may be tracing either dif- 
ferent physical properties or different regions of the outflow 
system. While SiO may reflect the shock strength of the jet 
itself, the other tracers may reflect physical conditions of the 
ambient medium such as the density of ambient cloud (Blake 
et al. 1995). 

4. SUMMARY 

The NGC 1333 IRAS 4A region was observed using the 
VLA in the 1.3 cm and 6.9 mm continuum, the NH3 (2, 2) 
and (3, 3) lines, and the SiO w = OJ=l— !>0 line, with an 
angular resolution of ~0''4, to image the circumstellar struc- 
tures and outflows of the protobinary system. The continuum 
emission mostly comes from dust in the circumstellar disks 
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and the protostellar envelopes. The high-resolution (0''3) im- 
ages of the NH3 lines show the detailed structures of the disks 
and outflows around the protostars. The molecular line images 
of relatively low resolutions (l"-2") were used to investigate 
the properties of the protostellar outflows. The main results are 
summarized as follows. 

1 . The continuum maps show two compact emission sources 
and a clumpy extended structure surrounding them. They may 
include the circumstellar disks, protostellar envelopes, and the 
common envelope of the binary system. The steep SED sug- 
gests that most of the 1.3 cm and 6.9 mm continuum emission 
comes from dust. The mass of the dense molecular gas in the 
binary system, imaged by the millimeter-continuum interfer- 
ometric observations, is ~^2.6 A/©. The IRAS 4A1 protostar- 
disk-envelope system is more massive than the A2 system by a 
factor of ^6. 

2. The NH3 line maps of the region immediately around 
IRAS 4A1 show at least two components. One is probably the 
circumstellar disk, and the other may be the blueshifted south- 
ern outflow. The blueshifted and redshifted emission sources of 
the central (disk) component was not spatially separated, sug- 
gesting that the disk may be rotating relatively slowly. Analysis 
of the rotation curve, as what was done with the A2 disk, was 
not possible with the angular resolution of this work, and sug- 
gests that the Al protostar is less massive than the A2 protostar 

3. Considering the properties of the disks and the lengths 
of the outflows, IRAS 4A1 may be younger than A2, in terms 
of collapse age. We suggests that the less massive member of a 
binary prestellar core system may start the protostellar collapse 
earlier than its companion. 

4. The molecular outflows are brighter in the NH3 (3, 3) line 
than in the (2, 2) line. The overall structure of the NH3 out- 
flows is similar to that of the SiO jets/outflows, though there 
are interesting differences in small-scale structures and inten- 
sity distributions. 

5. The 1" resolution map of the SiO line shows the highly 
turbulent structures of the deflected part of the IRAS 4A2 
northeastern jet. The SiO line is brightest at a peak located ^7" 
from the starting point of the deflected part, suggesting that 
the enhancement of SiO in the jet takes ^100 yr after the pro- 
posed jet-core collision. The NH3 line maps show an ambient- 
velocity structure near the impact point, which may be a part 
of the obstructing dense cloud core. 

6. Comparisons of the SiO, NH3, H2CO, and HCN maps of 
the IRAS 4A outflows show an interesting chemical dichotomy. 
The intensity distribution of SiO is anti-correlated with those of 
the others. This trend suggests that NH3, H2CO, and HCN may 
trace the same bulk of gas in the molecular outflows, while SiO 
traces a gas component more directly related to the primary jet. 
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